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Abstract — Optimization of structures has always been a 
subject of continuous interest in the field of structural 
engineering. The amount of research work and publications in 
this field show various mathematical approaches adopted to 
effectively use materials used for construction. A novel 
iterative Perfect Nodal Position Search Method (PNPSM) for 
the evolution of optimum design of trusses is presented . This 
method does the analysis, design, generates the optimum 
topology, and arrives at minimum sizes of members, subjected to 
given stress and movement constraints. Efficiency of a member 
is determined at every level of iteration, end nodes are moved to 
improve the efficiency of the member. A forceful refinement of 
length of each member and hence the movement of end nodes 
leads to a refined geometry of the truss . Three parameters are 
identified corresponding to member properties determining the 
extend of modification of length and repositioning of nodes. The 
effect of maximum percentage of modification of length on the 
converge to the minimum weight of truss is studied. 

Index Terms — Structural Optimization, Topology, 
Sizing ,Shape, Trusses, Iterative Process, nodal position . 


I. INTRODUCTION 

Two dimensional trusses are basic and commonly used form 
of construction . Evolutionary ideas in optimizing the size of 
members, shape and topology lead to preservation of 
precious material. 

Kulkarni et. al [1] , have presented a mutation-based 
real-coded genetic algorithm (MBRCGA) for sizing and 
layout optimization of planar and spatial truss structures is 
presented . The standard deviation of design variables has 
been used as a key factor in the adaptation of mutation 
operators. The reliability of the algorithm has been 
investigated in sizing and layout optimization, with both 
discrete and continuous design variables. A hybrid 
real-parameter genetic algorithm has been developed by 
Hwang et. al. [2] to solve optimization problems. The 
performance of the algorithm in discrete sizing variables and 
continuous configuration variables, both individually and 
combined has been studied. Tang et. al [3] have developed an 
Improved Genetic Algorithm for design optimization of truss 
structures with sizing, shape and topology variables with 
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mixed coding of integer and float types of variables. Sizing, 
geometry and topology optimization of trusses via force 
method and genetic algorithm has been introduced by 
Rahami et. al [4] which used a combination of energy and 
force methods in optimizing the weight of trusses.. Pavel et. 
al. [10] have presented a method for the simultaneous 
topology and size optimization of 2D and 3D trusses using 
evolutionary structural optimization with regard to commonly 
used topologies. An interesting Imperialistic Competitive 
Algorithm for truss structures has been presented by Hadi et. 
al. [8] inspired from social human phenomenon , in which 
some empires with lowest cost are considered the best and the 
rest of the countries in the neighborhood are considered 
colonies. The power of a country is inversely proportional to 
its cost. This has been used as a function to ultimately solve 
the optimization problem using Genetic Algorithm. 

II. OPTIMIZATION OF TRUSSES 

A Simultaneous Analysis and Design of Trusses 

Analysis of a truss involves, sequentially, the initialization 
of the problem with nodal positions, member sizing, 
connectivity, support conditions and external loads. In a 
Standard Direct Stiffness Matrix Method, a Global Stiffness 
Matrix is assembled from individual element stiffness 
matrices. Support conditions are imposed by manipulation of 
corresponding elements in global stiffness matrix. With the 
force vector on the right hand side, simultaneous equations 
are solved to get displacements at nodes and forces in 
members. Stresses in the material are determined to check the 
adequacy of cross sectional areas provided for all the 
members. 

B Sizing Optimization 

Sizing is mainly governed by the permissible stresses in the 
material under different conditions. The tensile stress limit 
could be a factor of the yield stress. But the slenderness ratio 
and Young’s modulus play a critical role in compressive 
stress limit. 

If , in a member, the stress in the material exceeds the 
permissible stress, the member is unsafe and hence the design 
is unacceptable. The cross sectional areas of such members 
need upward revision. At the same time, some members may 
have been oversized initially. The cross sectional area of 
those members need downward revision. This procedure is 
called sizing optimization. An ideal situation is the one for 
which the material strength gets fully exploited for all the 
members, simultaneously. 
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C Shape Optimization 

It involves finding out a shape of a truss which gives minimum 
weight. This is achieved by moving the nodes of a truss in an 
efficient manner such that the shape evolved is the best to 
effectively transfer the loads expected on it, at various nodes. 

D Topology Optimization. 

In trying to find out the perfect shape of truss, we may find 
that some members may turn out to be totally inefficient, with 
the cross sectional areas demanded by size optimization 
reducing to minimum. Judicious removal of such members 
from the configuration of the truss, give rise to the best 
topology. 

III. EVOLUTION OF PERFECT NODAL POSITION 
SEARCH METHOD 

Motivation to develop the proposed Perfect Nodal Position 
Search Method (PNPSM) comes from the necessity of an 
algorithm that uses the basic principles of structural 
mechanics effectively, efficiently and intelligently to reduce 
material consumption in trusses in a systematic and 
sequential manner. Though sizing optimization is member 
based, PNPSM uses nodes as points of attraction in shape 
and topology optimization. 

The proposed method uses the Direct Stiffness Matrix 
Method for the analysis of the truss, with the given size of 
members, shape and topology of the structure, as starting 
values of variables . As usual, while assembling the Global 
Stiffness Matrix, the resolved components of the element 
stiffness values are algebraically added at the positions 
representing the end nodes of a member. Initial values of 
member properties like A (Cross sectional area) and E 
(Young’s Modulus of the material) and 1 (Length of the 
member) are taken as inputs for the determination of element 
stiffness. 

Solution is obtained for forces and the actual axial stresses are 
determined. The ratio of actual stress to the permissible 
stress in the material at a cross section of a member is called 
utility ratio of that member , U. If U > 1 the member is 
unsafe. 

At a typical node, let us assume that there are 4 m’ number of 
members connected to it. Let the cross sectional areas be A 1? 
A 2 , A 3 ,.A m and lengths be -Im • 

Performing the analysis and getting a safe design, need not 
lead to an optimum design. Let us assume the utility ratios of 
the members be Ui,U 2 ,U 3 ,...U m , all of them being less than 
unity for a safe design. Utility ratios of members is the 
guiding premise for the formulation of Perfect Nodal 
Position Search Method (PNPSM) that lead to optimum 
design through a sequential achievement of the following 
criteria 


A. Criterion-1 : Sizing optimization is achieved by forcing 
change in cross sectional areas of members to reach values 
that give utility ratios close to unity. 

Achieving this condition for all the members , necessitates 
search for strategic node location, which has prompted the 
criterion given below. 

B Criterion-2 : If a node can be moved to a new position 
such that all the members connected to it reach the utility 
ratios very close to unity simultaneously, then that position is 
deemed to be the perfect position for the node and hence the 
PNPSM. 

These conditions, when satisfied for all the nodes will 
result in all the members of the truss reaching utility ratios 
nearing unity. 

The cross sectional areas are derived from the forces in 
members in sizing optimization. Some of the members may 
end up with negligible cross sectional areas en-route the 
adaptive search for perfect nodal positions. This possibility is 
well addressed in Criterion -3 in the following section. 






C. Criterion-3 : In the search process, If a member area 
demanded is miniscule (Negligible Cross sectional area) the 
indication is that such members are ineffective and the truss 
can perform without that member in question, as a part of the 
current configuration. Adaptive search may also encounter a 
situation where all the members at a node are ineffective. To 
tackle this situation, Criterion -4 is formulated. 
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D Criterion-4: If all the members at a node are ineffective, 
then the truss can survive without that node, indicating it can 
be collapsed to any node it is linked to. 

This gives us an opportunity to remove such nodes from the 
truss and end up with a better shape. 


While Criterion-1 concentrates on members, the others 
concentrate much on nodes. 



Fig. 3. Criterion for Removal a Member 

NEGLIGIBLE node 



Fig. 4. Criterion for Removal of a Node 


IV MOVEMENT OF NODES 

The multi objective optimization of trusses with same 
material can now be viewed as finding out the perfect nodal 
positions such that all the members of the truss, subjected to 
given set of constraints 

a) have non negligible cross sectional areas 

b) have utility ratios equal to unity (ideal situation) 

c) have a set of lengths such that total weight of the truss 
W, is minimum 

Where 

W =Unit wt X Sum of ( Ai X lj) — Eq. 1 

i = 1 to m, where Ai is cross sectional area and li is the length 
of of the i th member 

Movement of node changes the lengths of members connected 
to the node. The reverse is also true. The change in lengths of 
members connected to a node, moves the node. 

To achieve the minimum, the problem is approached from 
three angles, simultaneously, for a member. 


i. Reduce the length, if cross sectional area is to remain the 
same 

ii. Reduce the cross sectional areas if the lengths is to 

remain the same. 

iii. Increase the effectiveness 

With these in mind, the lengths of members connected at a 
node are changed based on their relative qualifications. At a 
node , where ‘m’ number of members are connected, for 
every member three factors are identified which are factors 
dependent on its length, cross sectional area and utility ratio. 
Factor for weight consideration (Ci) 

C n = Wi / sum of ( Wi) — Eq. 2 

i = 1 to m , where Wi is the weight of the i th member 

Factor for area consideration ( C 2 ) 

C 2i = Ai / sum of ( Ai) — Eq. 3 

i = 1 to m, where Ai is sectional area of the i th member 

Factor for inefficiency (C 3 ) 

C 3i = 1 - Ui ~~ Eq. 4 

where Ui is utility ratio of the i th member 

Total forced change in length of ith member ,dl in the 
direction of member 

dli = MF . Cii . C 2i . C 3i . li —- Eq. 5 

where MF is the desired maximum percentage modification 
desired per iteration 



Fig. 5. Forced Change in Length of a Member 

When the member length changes by dl , the node will be 
moved with respect to its original position by dxl in the global 
X axis and dyl in the Global Y axis. 

X new = X old + sum of (dxl) for all members -- Eq. 6 

Y new = Y old + sum of (dyl) for all member — Eq. 7 

The factor C2, defined for the area of a member becomes 
negligible if the member is ineffective. The factor C 3 which 
is meant for the inefficiency of a member reduces to zero, 
when utility ratio is unity. This means , the length of a member 
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is not altered in both the cases. If all the members at a joint are 
efficient, the node is not moved. This ensures the convergence 
of the solution for optimization. On the other hand, the node 
will be forced to move relatively when every member is 
subjected to a forced change in length, every time checking 
the node position for a possible set of movement restrictions 
specified in the problem. This is repeated for all the nodes and 
the changed configuration of the truss is recorded. 

V ITERATIVE PROCEDURE 


P = 10 kips ( 44.48 kN) 
a = 120 in (3048 mm) 

density = 0.1 lb/ cu. in ( 2768 kg/ Cu. m) 
e = 10000 ksi ( 68947.6 N/ Sq. mm ) 
limit = 25 ksi ( 172.0 N / Sq. mm ) 


Table 1. Constraints for Nodal Movements 


The PNPSM is an iterative procedure with distinct loops for 
sizing and shape optimization. Topology optimization is 
achieved during the course of shape optimization. 

Starting with the initial geometry, member properties, 
loading and support conditions and constraints, solution is 
obtained for the stresses. Cross-sectional areas are increased 
or decreased iteratively to obtain a safe sizing optimization 
for the shape and topology, till the utility ratios stabilized. 
This is named as the sizing loop. 

Member lengths are modified as per Eq. 5 to effect change 
in nodal positions and the whole procedure is repeated to get 
another set of stabilized utility ratios. This is named as the 
combined shape optimization loop. Size optimization loop is 
a part of shape optimization loop. The procedure is repeated 
till we get a stabilized set of utility ratios of all effective 
members equal to unity. 

While iterative loop in progress, some of the members of the 
truss are identified ineffective and Young’s modulus values 
of such members are considered negligible for the 
consecutive loop. The final shape of the truss without these 
ineffective members is the optimum topology. 

Fig. 7 shows the flowchart for the implementation of the 
PNPSM. 


JOINT MOVEMENT C ONSTRAINTS 

JOINT No. 

COORDINATES 

Permissible Freedom 


X 

Y 

MinX 

MaxX 

MinY 

MaxY 

1 

0 

120 

0 

0 

120 

120 

2 

120 

120 

100 

140 

100 

140 

3 

240 

120 

220 

260 

100 

140 

4 

360 

120 

360 

360 

50 

90 

5 

0 

0 

0 

0 

0 

0 

6 

120 

0 

100 

140 

-20 

20 

7 

240 

0 

220 

260 

-20 

20 

8 

360 

0 

360 

360 

20 

60 


Additional Conditions : X 6 = X 2 , X 7 =X 3 , X 8 = X 4 = 360 
Fig. 6. Benchmark Problem 


VII RESULTS AND DISCUSSION 

The optimum topology evolved is shown in Fig. 8. Table 2 
shows the set of Length, Area and Utility ratio for every 
member at the instance of optimum design. It is noted that 
PNPSM clearly identifies member numbers 3,7,8,9 and 15 as 
ineffective. The positional changes of the nodes for the 
optimum configuration have been affected only on the 
foundation of utility ratio wherein permissible stresses both 
in tension and compression remain the same, as stated in the 
benchmark problem. 


VI NUMERICAL EXAMPLE 

The 15 bar truss problem, shown in Fig.6, solved by many 
researchers [1] has been treated as benchmark to check the 
efficiency of the PNPSM algorithm. The optimum design is to 
be achieved with the properties and movement restrictions 
stated in Table 1. 



effective member 
ineffective member 



p 


Fig 8. Configuration for Minimum Weight - 
Topology 

Table 3 shows the sizing and layout variables obtained by 
PNPSM in comparison with the results given in references. 
It is seen that the results obtained are in agreement and 
showing a further improvement in optimum design 

Fig. 9 shows the weight reduction of the truss corresponding 
to iterations performed with a typical Modification Factor 5% 

A study has been conducted to know the effect of 
Modification Factor (MF) in PNPSM on the convergence to 
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the optimum design of the 15-Bar Truss. Fig. 10 shows the 
individual convergences for various Modification Factors. It 
is seen that the value of Modification Factor up to 15% has no 
much effect on the convergence to the optimum design in 
terms of the shape, topology, sizing and finally, the weight. 
Fig. 11 shows the effect of Modification Factor on the 
Convergence to minimum weight. 


Perfect Nodal Position Search Method as presented and 
demonstrated clearly has an edge over other techniques 
discussed . The method is appealing as convergence is faster 
and the search for optimum shape, topology and member 
sizing is accomplished simultaneously. Suitable 
modifications to element stiffness matrix, extends its usage 
potential to 3D trusses too. 


VIII CONCLUSIONS 





Fig 8. Configuration for Minimum Weight - Topology 


Fig. 7. Flow Chart for Perfect Nodal Position Search Method 
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Table 2. Final Cross Sectional Areas and Utility Ratios 


RESULTS WITH REDUCTION FAC TOR = 5 % 

15 MEMBER TRUSS OPTIMUM DESIGN FORC ES ON MEMBERS 


Length 

AREA 

WT 

FORCE 

Type 

STRE 

SS 

Allow. 

Stress 

RATIO 

Memb. No. 

in 

Sq. in 

kips 

kips 


ksi 

ksi 


1 

120.0000 

0.8671 

0.0104 

22.060 

Tens 

25.00 

25.00 

1.00 

2 

119.9030 

0.7321 

0.0088 

17.750 

Tens 

25.00 

25.00 

1.00 

3 

123.8110 

0.0000 

0.0000 

0.000 

INEFFECT 

0.00 

25.00 

0.00 

4 

120.0000 

1.1332 

0.0136 

-27.940 

Comp 

-25.00 

25.00 

1.00 

5 

119.9030 

0.4670 

0.0056 

-11.910 

Comp 

-25.00 

25.00 

1.00 

6 

121.7860 

0.4050 

0.0049 

-9.840 

Comp 

-25.00 

25.00 

1.00 

7 

120.0000 

0.0000 

0.0000 

0.000 

INEFFECT 

0.00 

25.00 

0.00 

8 

120.3080 

0.0000 

0.0000 

0.000 

INEFFECT 

0.00 

25.00 

0.00 

9 

70.0000 

0.0000 

0.0000 

0.000 

INEFFECT 

0.00 

25.00 

0.00 

10 

169.7060 

0.4700 

0.0080 

11.230 

Tens 

25.00 

25.00 

1.00 

11 

169.7060 

0.0951 

0.0016 

-2.910 

Comp 

-25.00 

25.00 

1.00 

12 

169.7860 

0.0960 

0.0016 

3.280 

Tens 

25.00 

25.00 

1.00 

13 

169.7060 

0.4710 

0.0080 

-11.660 

Comp 

-25.00 

25.00 

1.00 

14 

156.3420 

0.5200 

0.0081 

12.580 

Tens 

25.00 

25.00 

1.00 

15 

150.2050 

0.0000 

0.0000 

0.000 

INEFFECT 

0.00 

25.00 

0.00 

TOTAL WEIGHT OF TRUSS 

0.07066 

kips 



Table 3. Comparison of Results for Benchmark Problem. 


15- MEMBER TRUSS - OPTIMUM DESIGN - C OMPARISON OF RESULTS 


EARLIER WORKS [1] 

Present 

Work 

MEMBER 

Gholizadeh 

Tang et 

Hwang and 

Rahami et 

Kulkarni et 

} AREA l 

et. al. [12] 

al. [3] 

He [2] 

al. [41 

al. [1] 

Sizing variables (in. 2 ) 

A1 

0.954 

1.081 

0.954 

1.081 

0.954 

0.8671 

A2 

0.539 

0.539 

1.081 

0.539 

0.539 

0.7321 

A3 

0.27 

0.287 

0.440 

0.287 

0.111 

0.0000 

A4 

1.081 

0.954 

1.174 

0.954 

0.954 

1.1332 

A5 

0.539 

0.954 

1.488 

0.539 

0.539 

0.4670 

A6 

0.174 

0.220 

0.027 

0.141 

0.347 

0.4050 

A7 

0.111 

0.111 

0.270 

0.111 

0.111 

0.0000 

A8 

0.111 

0.111 

0.347 

0.111 

0.111 

0.0000 

A9 

0.44 

0.287 

0.220 

0.539 

0.111 

0.0000 

A10 

0.44 

0.220 

0.440 

0.440 

0.440 

0.4700 

All 

0.347 

0.440 

0.347 

0.539 

0.44 

0.0951 

A12 

0.22 

0.440 

0.220 

0.270 

0.174 

0.0960 

A13 

0.22 

0.111 

0.270 

0.220 

0.174 

0.4710 

A14 

0.174 

0.220 

0.440 

0.141 

0.347 

0.5200 

A15 

0.27 

0.347 

0.220 

0.287 

0.111 

0.0000 

Layout variables (in.) 

X2 

113.65 

133.612 

118.346 

101.5775 

105.7835 

120.000 

X3 

254.47 

234.752 

225.209 

227.9112 

258.5965 

239.903 

Y2 

128.97 

100.449 

119.046 

134.7986 

133.6284 

120.000 

Y3 

115.73 

104.738 

105.086 

128.2206 

105.0023 

120.098 

Y4 

59.364 

73.762 

63.375 

54.8630 

54.4546 

90.000 

Y6 

-12.733 

-10.067 

-20.000 

-16.4484 

-19.9290 

0.000 

Y7 

3.5467 

-1.339 

-20.000 

-13.3007 

3.6223 

-0.211 

Y8 

59.29 

50.402 

57.722 

54.8572 

54.4474 

20.000 

Weight (lbs) 

73.93 

79.820 

104.573 

76.6854 

72.5152 

70.660 
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Fig 9. Typical Graph Showing Weight reduction of Truss 



Fig 10. Effect of Modification Factor on The Optimum Design of the 15-Bar Truss 




*MF= i % 
-*-\IF= 2% 
-*-\IF= 5% 
-#-\IF= 10% 
-*-\IF= 15% 


ITERATIONS 


Fig 11. Effect of Modification Factor on Convergence to minimum weight 
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